Interplay between lattice and spin states degree of freedom in the FeSe 
superconductor: dynamic spin state instabilities 
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Polarized Raman-scattering spectra of superconducting, single-crystalline FeSe evidence pro- 
nounced phonon anomalies with temperature reduction. A large (~ 6.5%) hardening of the Bi g (Pe) 
phonon mode is attributed to the suppression of local fluctuations of the iron spin state with the 
gradual decrease of the iron paramagnetic moment. The ab-initio lattice dynamic calculations sup- 
port this conclusion. The enhancement of the low-frequency spectral weight above the structural 
phase transition temperature T 3 and its change below T s is discussed in relation with the opening 
of an energy gap between low (S = 0) and higher spin states which prevents magnetic order in 
FeSe. The very narrow phonon line widths compared to observations in FeTe suggests the absence 
of intermediate spin states in the fluctuating spin state manifold in FeSe. 

PACS numbers: 74.70.Xa, 74.25.Kc 



I. INTRODUCTION 

The iron selenide FeSe superconductor, being simplest 
with respect to its crystal structure, is the most unusual 
compound among the large family of iron based high tem- 
perature superconductors (HTSC). In fact, it is the only 
compound which remains in the nonmagnetic metallic 
state for temperatures down to 0.02 K»i Furthermore, 
it becomes superconducting in the pure stoichiometry 
phase^i^ i.e. without doping in contrast to other HTSC 
members. At first glance it is a conventional supercon- 
ductor with a modest transition temperature of T c ~ 8 
K at ambient pressured However, it shows an enormous 
growth of T c up to 37 K at pressures of 9 GPa^— sug- 
gesting that the lattice subsystem plays a role in the su- 
perconductivity. 

From magnetic point of view FeSe and its main 
building blocks of FeSe4 tetrahedra demonstrate amaz- 
ing features: A stoichiometric FeSe thin film on a c- 
plane sapphire substrate is a ferromagnetic metal with 
a Curie temperature above 300 K4 Moreover, in the 
recently discovered iron selenides intercalated by alkali 
ions Mo.8+ !E Fei.6+„Se2 (M = K, Rb, CsjSriS the FeSe 4 
tetrahedra in the FeSe layers stay in two different spin 
states. In the insulating and magnetically ordered part 
of the phase separated sample the FeSe4 tetrahedra are 
in the high spin state with a magnetic moment of 3.3 (j,b 
per tetrahedron at 10 KM- In the nonmagnetic metal- 
lic and superconductive phase of the Mo_s+xF^i.6+y^2 
system the FeSe4 tetrahedra remain in their low spin 
state (S = 0). All of these observations point to an 
easy switching between the low (S = 0), intermediate 
(5 = 1), and high (S = 2) iron spin states in the layers 
of FeSe4 tetrahedra and a pronounced instability of the 



iron spin states. This instability is governed by distor- 
tions of the tetrahedra and in particular by the flexible 
change of the height (z-coordinate) of the Se ions above 
the Fe plane. Thus, it suggests the presence of rather 
specific longitudinal spin fluctuations in FeSe which are 
intimately connected to the lattice. 

Actually, nonuniform, 5 / 0. antiferromagnetic spin 
fluctuations in FeSe, seen in NMR experiments^ are re- 
ported to be strong and increasing near the supercon- 
ducting transition temperature. Furthermore, applied 
pressure enhances these spin fluctuations together with 
an increasing T C M However, the uniform q = paramag- 
netic spin susceptibility, which is derived from the NMR 
Knight shift, shows an almost linear decrease under tem- 
perature lowering^ This observation can be interpreted 
as an unexpected suppression of the uniform spin fluctu- 
ations in FeSe with decreasing temperature. 

The temperature evolution of spin fluctuations should 
be seen in the lattice dynamics through a high sensi- 
tivity of the phonon energies on the varying iron spin 
states^— Already earlier first principle studies for 
"1111" and "122" systems have revealed a strong depen- 
dence of the z-atom positions from the Fe spin stated It 
was shown that under external pressure, the high Fe-spin 
state structure should transform to a structure with low 
Fe-spin state and significantly reduced c-axis. Recent 
experimental and theoretical studies of the zone center 
phonons in the isostructural FeTe - a magnetic and non- 
superconductive relative of FeSe - have shown that the 
real magnetic moment of the Fe ions should be included 
in spin polarised lattice dynamic calculations to achieve 
an agreement between theory and experiment i 19 ' 20 Fur- 
thermore, the high frequency part of phonon spectra 
which in the iron chalcogenides HTSC mostly consists of 
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iron vibrations shows a giant high frequency shift under 
a decrease of the iron magnetic moment. The underlying 
physics of this phenomenon is related to strong orbital 
reordering under a change of the spin state. 

Here we report results from Raman light scattering and 
ab-initio lattice dynamic calculations of recently grown 
high quality FeSe single crystals. We demonstrate that a 
dynamical spin state instability is realized in FeSe at high 
temperatures. Strong dynamical fluctuations of the iron 
spin state are seen as a renormalization of the B lg (Fe) 
phonon energy at high temperatures with a dominance 
of the high spin state. Under temperature lowering these 
fluctuations decrease due to the appearance of an energy 
gap between different spin states with mainly low spin 
state presence. This is observed as an additional phonon 
frequency hardening which highly exceeds the usual tem- 
perature dependent hardening given by the Griineisen 
law. The presence of mainly low and high spin fluctuating 
states in the paramagnetic ground state in FeSe is seen 
as a one order narrowing of the phonon linewidth com- 
pared to the width of those modes in FeTe where a highly 
orbital frustrated intermediate spin state determines its 
properties^ The suppression of spin state fluctuations is 
continuous and smooth with lowering temperatures. Be- 
low 150 K we observe a highly polarized low frequency 
background which becomes a gap-like structure. This 
feature is attributed to a gap opening between a manifold 
of intermediate spin states and the paramagnetic ground 
state in FeSe. Our findings shed light on the mechanism 
that prevents magnetic order and the appearance of a 
superconducting state in FeSe. 

II. EXPERIMENTAL DETAILS 

We have studied FeSe single crystals grown in evac- 
uated quartz ampoules using a KCI/AICI3 fluxi^i The 
very well characterized samples were from the same se- 
ries as in previous investigations ! 21 ^ 22 The present sam- 
ples composition is FeSeo.96- Raman spectra were mea- 
sured on freshly cleaved surfaces in a quasi-back scatter- 
ing configuration using a A = 532.1 nm solid-state laser. 
The laser power was set to 5 mW with a spot diameter 
of approximately 100 /jm to avoid heating effects. All 
measurements were carried out in a closed-cycle cryostat 
(Oxford/Cryomech Optistat) in the temperature range 
from 7 to 295 K. The spectra were collected via a triple 



spectrometer (Dilor-XY-500) by a liquid nitrogen cooled 
CCD (Horiba Jobin Yvon, Spectrum One CCD-3000V). 
The polarization configuration of the Raman spectra is 
denoted by (e,e s ). We used a- and o-notation for the 
tetragonal crystallographic axes and x- and y-notations 
for the orthorhombic axes which are rotated by 45° from 
the a- and 6-axes. 

III. RESULTS 

At room temperature, FeSe has the tetragonal PbO 
structure (space group P4/nmm)^ In this phase Fe and 
Se ions occupy 2a and 2c Wyckoff positions, respectively. 
Symmetry analysis shows that there are four Raman- 
active modes [r flama „ = Ai g (Se) + Bi g (Fe) + 2E g (Se, 
Fe)] and two infrared-active modes [Tjh = A 2u (Se, Fe) 
+ E„(Se, Fe)]. The Raman tensors are given by: 

fa 0\ fc 0\ / -e\ 

A lg = a , B lg = -c , and E g = e . 
\0 bj \0 0/ \-e e J 

Indeed, two sharp peaks are observed in the tetragonal 
phase at room temperature at 179.8 cm -1 [Ai g (Se)] and 
193.9 cm -1 [Bi g (Fe)], see the (aa) spectrum in Fig. Q] 
On the one hand, the narrow linewidth of the peaks indi- 
cates a high quality of the samples. On the other hand, 
a monotonic phonon mode sharpening was observed in 
Fei+yTei-^Sea; with increasing Se content^ This fact 
is quite surprising since one may expect a concentra- 
tion dependent broadening of the phonon lineshape due 
to the Se substitution-induced disorder. It was argued 
that spin-phonon coupling may be responsible for this 
behaviour. 

At T s = 90 K, FeSe undergoes a second-order fer- 
roelastic phase transition from the tetragonal to an or- 
thorhombic phase, and below 90 K, its space symmetry 
group is Cmma, with a y/2x v2 supercell enlargement in 
the basal planed In the low-temperature orthorhombic 
phase Fe and Se ions occupy 4a and 4o Wyckoff positions, 
respectively. Symmetry analysis gives six Raman-active 
modes [r ' Raman = A g (Se) + B lg (Fe) + 2B 2g (Se,Fe) + 
2B 3g (Se,Fe)] and six IR-active modes [T IR = 2Bi„(Se,Fe) 
+ 2B 2u (Se,Fe) + 2B 3u (Se,Fe)]. The corresponding Ra- 
man tensors are given by: 



fa 0\ fO d 0\ fO e\ /0 0\ 

A g = b , B lg = d , B 2g = Land B 3g = / . 

yooc/ \000y \eooy \o / oy 



Figure [2] presents polarized Raman spectra of FeSe at 
T = 7 K. The Ai g mode in the tetragonal structure 
changes into an A g mode in the orthorhombic structure 
by the compatibility relation of the group theory and is 



observable at 182.6 cm -1 in (aa) and (xx) scattering con- 
figurations. The Bi g mode at 206.9 cm -1 is observed in 
(aa) and (xy) scattering geometry. Besides, weak lines 
with frequencies around 236, 290 and 357 cm -1 appear 
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FIG. 1. (Color online) Polarized Raman spectra of FeSe at 
room temperature. The inset shows the atomic displacement 
patterns of the Raman-active phonon modes. 



FIG. 2. (Color online) Polarized Raman spectra of FeSe at 
T = 7 K. The left inset shows a comparison between FeSe 
and FeTe at T = 50 K. The right inset shows the temperature 
evolution of the (aa) polarized spectra. 



in (xx) polarized spectra at low temperatures, which can 
be assigned to B2 ff /B3 9 modes. The frequencies of these 
modes have no counterparts in Raman active phonon 
modes of FeSe, therefore we assign them to a small ad- 
mixture of iron/selenium oxides £ Here we want to high- 
light that the positions of the phonon lines in our inves- 
tigations do not coincide with those of a Raman study in 
Refi2£. The reason may be related to sample quality of 
the previous investigation. A much better agreement is 
found with the data of a phonon density of states study 
by 57 Fe nuclear inelastic scattering^ Here the B l£ ,(Fe) 
phonon mode was identified at the frequency of 196.8 
cm -1 . 

The temperature evolution of the (aa) phonon spectra 
is presented in the right inset of Fig. [5J To get more in- 
sight into the phonon dynamics, we analyzed the temper- 
ature dependence of the phonon lines parameters. Lines 
were fitted with Lorentzians in the whole investigated 
temperature range and results are presented in Fig. [3] 
where the temperature dependent parameters of the A lg 
and Bi 9 phonon modes are plotted. 

The T dependencies of the frequencies of both modes 
are smooth, see Fig. [3J (a) and (b), their positions shift 
to a higher energy with the temperature decreasing. The 
Ai s (Se) mode hardens by 3 cm -1 while the B lg (Fe) mode 
undergoes a large hardening by 14 cm -1 (6.5%) in the 
whole temperature interval. An analysis shows that such 
a behaviour can not be explained only by the crystal 
lattice contraction upon cooling. Using the data for the 
lattice parameters from Ref3£, we can estimate that in 
the frame of the Gruneisen law, the parameter 7, must 
be equal to ~ 1.04 for the Ai ff (Se) mode and ~ 5.09 
for the B lff (Fe) mode. Both of these values deviate from 
the conventional value 7 ~ 2 as well as from 7Ai g (Te) = 
2.85 and "fBig(Fe) = 2.15 as observed by us in FeTeJ^ A 
hardening of the phonon density of state in the region of 
the zone centred Bi 9 (Fe) phonon energy was observed in 
Refi2£ and can be estimated with a Gruneisen parameter 



Of iBlg(Fe) ~ 4. 

The first value is quite reasonable. However, the latter 
deviates from the conventional values. The hardening of 
all phonon modes under pressure was observed in Refill 
and the estimated Gruneisen parameter is jBig(Fe) = 
0.9. Given that only a part of the observed frequency 
shift can be attributed to 1% of the unit cell volume ther- 
mal contraction 2 ^ and to an increasing anharmonicity 
of the vibrational modes with increasing temperature^ 
we believe that the anomalous Bi g (Fe) phonon harden- 
ing is connected with the spin-state of Fe ions changing 
as well. 

The solid lines in the Fig. [3J (c) and (d) are a fit of the 
temperature dependence of the phonon linewidths using 

r(r) = r fi+ % ) (l) 

with To = 1.52/1.98 cm -1 , the eigenfrequency ujq = 
182.6/206.9 cm -1 , and dj being a mode dependent fit 
parameter for Ai s /Bi g phonon modes. The data can 
be described very well down to the lowest temperature 
by this expression, which includes cubic anharmonicity 
effects . 31 i 32 Noteworthy, however, is the superimposed 
change of curvature leading to an additional drop of 
linewidth for temperatures around T s for the Bi s phonon 
mode. 

Figures [3] (e, f) present the integrated intensity of two 
phonon lines, in analogy to the frequency and linewidth 
of the phonon modes. It is noteworthy that the inten- 
sity evolution is different for each mode. However, both 
curves have features which are consistent with the tem- 
perature of the structural phase transition, T s = 90 K, 
and the temperature T* = 230 K, at which a new fea- 
ture, associated with changes in the structural sector, 
has been found in studies of Refj^ 2 .. Evidently also, the 
Ai g (Se) mode reveals a regular deviation from a purely 
Gruneisen behaviour for T < 90 K (Fig. [3J (a)). 
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FIG. 3. (Color online) Temperature dependence of the Ai 9 (Se) and Bi 9 (Fe) phonon line parameters in FeSe. (a, b) - phonon 
frequencies. Black dashed lines correspond to Griineisen dependencies with 7ai 9 (s e ) = 1-04, 7si 9 (Fe) = 5.09, and 7si 9 (Fe) = 0.9. 
(c, d) - linewidths (FWHM). Black dashed lines correspond to fits with Eq. (1). (e, f) - integrated intensities. 



As mentioned above, enhanced spin fluctuations at the 
nesting vector may cause a spin-density wave. A previous 
Raman scattering study of BaFe2As2^ revealed changes 
of the electronic properties at the SDW transition. The 
nodal and anti-nodal gap excitations with different sym- 
metry were observed in the spectra and theoretically an- 
alyzed. In order to clarify the trend in going from Fe-As 
based materials to FeSe, we performed a detailed study 
of polarization and temperature dependence of our spec- 
tra. The temperature dependence of the Raman spectra 
in (aa), (ab), (xx), and (xy) polarization configurations 
is shown in Fig. @] 

At all temperatures the spectra have a low-energy 
quasi-elastic spectral weight except in (ab) scattering ge- 
ometry, where the quasi-elastic signal appears only at 
temperatures below ~ 150 K. However, the temperature 
dependence of the low-frequency part is very different 
for (aa), (ab), (xx), and (xy) polarizations. In (aa) scat- 
tering geometry, the quasi-elastic signal only slightly de- 
creases with temperature reduction, in (xy) it is nearly 
constant, and in (xx) it has a slightly different shape 
at low temperatures. Drastic changes are observed in 
the (ab) geometry. The low-energy scattering intensity 
abruptly increases at temperatures below 150 K and then 
a gap-like structure forms together with a very broad 
hump at 350 cm -1 below T s . Simultaneously, the low- 
energy quasi-elastic intensity starts to evolve. Beside, 
a line with the frequency 182 cm' 1 appears in the (ab) 
spectra below T s due to orthorhombic distortions and the 
appearance of twins as well as small twins rotations in 
a6-plane. Figure [5] (a) shows difference spectra between 



various temperatures and 150 K. The zero levels are in- 
dicated by the horizontal lines with the same color as the 
spectra. Fig. [5] (b) demonstrates a temperature depen- 
dence of the integrated intensity of the quasi-elastic tail 
which has a characteristic maximum at T s , and is there- 
fore related to the second order structural phase transi- 
tion. However, it does not show a tendency to go to zero 
at low temperatures. Fig. [5] (c) shows the temperature 
evolution of the frequency of the most intensive part of 
the hump. 



IV. DISCUSSION 

At first, two issues of our observations should be ad- 
dressed: (i) the strong hardening of the Bi 9 (Fe) phonon 
frequency which is more than two times stronger com- 
pares to the expected one from lattice shrinking and the 
Griineisen law and also very different from the modest 
hardening of the Ai g (Se) phonon mode; and (ii) the ex- 
istence of a very narrow phonon linewidth on the back- 
ground of the strong and continuous phonon hardening, 
meaning that the mechanism of hardening does not con- 
tribute to the phonon decay. Note, that in the FeTe we 
observed a very conventional hardening of the Bi g (Fe) 
mode with 7sig(Fe) = 2.15 which is close to the expected 
7 ~ 2J£ Hence, there is no specific mechanism of hard- 
ening in the FeTe. At the same time both Raman active 
modes in FeTe demonstrate a very unconventional line 
width, one order larger compared to FeSe. 

The difference between these two compounds is the 
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FIG. 4. (Color online) Temperature dependent Raman spec- 
tra of FeSe. The low energy quasielastic spectral weight does 
not show a pronounced temperature evolution in the (aa), 
(xx) and (xy) channels of scattering. In the (ab) channel 
(B2 9 channel in the tetragonal setting) quasielastic scattering 
is not present at high temperature but develops below 150 K. 



spin state of iron. In Fe 1+J/ Te for temperatures be- 
low magnetic ordering the magnetic moment is around 
2 (Us/Fe depending on the excess of iron content 
This magnetic moment supposes the iron spin state to 
be intermediate 5 = 1 which is highly orbital frustrated. 
Indeed, a few types of alternating orbital ordering pat- 
terns, which are consistent with the magnetically ordered 
structure, have been found to occur in iron pnictide su- 
perconductors^ In our paper— we argued that in FeTe 
the orbital frustration mechanism resulting from the Fe 
intermediate spin state creates an additional channel for 
phonon decay. As a proof of this we observed an al- 
most doubling of the Bi 9 (Fe) phonon linewidth at the 
structural/magnetic phase transition at which strong or- 
bital fluctuations are naturally expected. In contrast, in 
FeSe no visible response of linewidth of both B lfl (Fe) and 
Ai 9 (Se) phonon modes at the structural phase transition 
is observed. Thus, one can safely conclude that the or- 
bital channel of phonon decay is not dominant in FeSe 
and, particularly, it is not relevant to the mechanism of 
the structural phase transitions in FeSe. Note, that at 
the PA/nmm — > Cmma second order ferroelastic phase 
transition, which is realized in FeSe, the optic phonons 
are not affected and some acoustic mode demonstrates 
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FIG. 5. (Color online) (a) Differential Raman spectra from 
150 K in (ab) polarization. The zero levels are shown by 
the horizontal lines with the same color as the spectra, (b) 
Integrated intensity of the quasi-elastic tail, (c) Energy of the 
hump as function of temperature. 



a soft mode behaviour through lowering of its sound ve- 
locity for wave vectors much lower than the Brillouin 
zone boundary. In FeSe this effect has been studiedi 22 i 36 
Such kind of acoustic phonon softening also does not af- 
fect the standard lattice anharmonicity channel of optic 
phonon decay. Therefore the phonons linewidths remain 
unchanged through passing T s (see Fig.[3]c, d). The only 
manifestation of the structural phase transition in FeSe 
in our Raman spectra is seen in the continuous increase 
in scattering intensity below T s (see Fig. [3] e, f). This 
effect can be connected with some rearrangement in po- 
larizability caused by orthorhombic distortions and its 
increase under lowering temperature. 

Concluding the experimental observations, the above 
comparison of FeTe and FeSe phonon is compatible with 
an intermediate spin state for FeTe in the whole investi- 
gated temperature regime. As this state is unaltered it 
does not contribute to phonon hardening, while it is a 
source of strong orbital fluctuations. On the other hand, 
the mechanism of strong phonon hardening for FeSe does 
not contain frustrating orbital degrees of freedom. At low 
temperature, however, the iron ions should be in the low 
spin state which has no spin orbital frustration similar as 
the high spin states. We naturally arrive to the conclu- 
sion that at high temperatures the iron ions in FeSe reside 
in the high spin state and the source of the strong B lg (Fe) 
phonon hardening is a dynamical crossover from the high 
to low spin state with temperature lowering. The corre- 
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sponding strong hardening of the high frequency part of 
the phonon spectrum in FeTe upon decrease of iron mag- 
netic moment has already previously been found in our 
earlier studies^ 

To further illustrate this interplay we performed ab- 
inito lattice dynamic calculations of FeSe under varying 
magnetic moments at the iron sites within the framework 
of density functional theory (DFT) and a frozen phonon 
approach. We applied the all-electron full-potential 
linearized augmented plane wave method (Elk code)^l 
with the local spin density approximation (LSDA}2£ for 
the exchange correlation potential and with the revised 
generalized gradient approximation of Perdew-Burke- 
Ernzerhof (PBEsol)^ The calculations were performed 
on a 9 x 9 x 6 grid which corresponds to 60 points in 
the irreducible Brillouin zone. We used the experimen- 
tal unit cell parameters for FeSe in the tetragonal phase 
at 300 K£ at ambient pressure with the unit cell lattice 
constants fixed. To reproduce the fluctuating magnetic 
moments in the paramagnetic phase we replaced them by 
static ones and suggested some kind of magnetic order. 
The spin-polarized calculations were performed with and 
without structure optimization for M^e = and with 
structure optimization for Mjr e = 2/is- In the latter case 
a checkerboard type of antiferromagnetic ordering with 
magnetic moments along the c-axis was supposed. The 
results of our calculations and a comparison with the ex- 
periment are given in Table I. 

As we expected the results for optimized structures do 
not coincide with the experiments giving highly overesti- 
mated Aip(Se) and Bi s (Fe) phonon frequencies for both 
M = and M = 2 fig iron magnetic moments. This is due 
to a too low value of the calculated Se ion height (zse) 
that correspondingly decreases the Fe-Se distance com- 
pared to the experimental ones. However, these results 
simulate the phonon frequency behavior under decrease 
of iron magnetic moments. As one can see, the frequency 
of the Bx s (Fe) phonon drastically increases when the size 
of magnetic moments changes from M = 2/xstoM = 0, 
while the Ai ff (Se) phonon frequency undergoes modest 
change. Note that we do not discuss here the question 
about the importance and necessity of spin-polarized cal- 
culations in iron HTSC to obtain true phonon frequen- 
cies and electron-phonon coupling constants that are now 
generally acceptedii but we focus on the phonon energy 
dependence on the size of the magnetic moment. 

The spin polarized calculations without structure op- 
timization with M = give much better coincidence with 
the experiment for the Ai a (Se) mode but again overesti- 
mate the Bi 9 (Fe) mode. This overestimation is expected, 
as the experimental value of z$ e at 300 K is higher than 
the calculated zs e not only for M = but also for static 
magnetic moment Mj? e = 2fj,s- The high value of zs e cor- 
responds to a high magnitude of the magnetic moment 
(see for instance^). Thus, one can conclude that under 
dynamical spin state fluctuations the iron paramagnetic 
moment at 300 K should exceed 2 This estimation is 
in qualitative agreement with probing of the local mag- 



netic moment in the FeTeo.3Seo.7 by Fe Kp x-ray emis- 
sion spectroscopy— At room temperature they found 2 
fiB for a fluctuating local magnetic moment in the para- 
magnetic phase. Note, that in the Raman studies of the 
superconductive Fei+yTei-^Se^ compounds the Bi g (Fe) 
mode shows almost the same features of frequency hard- 
ening and linewidth narrowing as in our case^ This im- 
plies the same mechanism of this phenomena for these 
compositions. 

Summarizing our treatment of experimental and cal- 
culated data on the width and anomalous hardening of 
the Bi ff (Fe) phonon we conclude that the paramagnetic 
ground state in FeSe can be presented as a superposition 
of low and high iron spin states with larger contribu- 
tion of the high spin state at room temperature. Un- 
der temperature lowering the comparative contribution 
of those states undergoes a rearrangement with a contin- 
uous growth of the low spin state part over the high spin 
state part. This is seen as a decrease of the ground state 
paramagnetic moment. A specific spin-state lattice cou- 
pling is contributing to the Bi ff (Fe) phonon hardening 
which does not create an additional channel for phonon 
decay and linewidth increase. Our indirect observation 
of the Fe paramagnetic moment decrease is in nice agree- 
ment with the almost linear lowering of the paramagnetic 
spin susceptibility with decreasing temperatures in FeSe 
which has been observed in the NMR experiments^ The 
idea of representing the magnetic ground state in iron 
based HTSC as a coherent superposition of different spin 
states has recently been proposed and discussed.— Here 
we conclude that the paramagnetic ground state in FeSe 
combines both low and high spin states and does not 
include the intermediate spin state. 

Furthermore, the result of our calculation for non op- 
timized structure with M = at 300 K can be used to 
roughly estimate the expected phonon frequencies at low 
temperatures. Indeed, in first approximation one can 
neglect contributions to the hardening from thermal lat- 
tice contraction as the main contribution comes from the 
change of the iron spin state. Then supposing M = 0, 
as it should be expected at lowest T, we still have a 
disagreement between theory (218.8 cm -1 ) and experi- 
ment (206.5 cm -1 ) in the Bi 9 (Fe) phonon frequency at 
7 K, supposing that M ^ 0. This is evidence for resid- 
ual spin fluctuations (nonzero contribution of high spin 
state) which survive even at 7 K - the lowest temper- 
ature of our experiment. In other words, FeSe remains 
in a paramagnetic metallic state upon approaching the 
superconductive state. Interestingly this estimation has 
resemblance with results of an inelastic neutron scatter- 
ing experiment in FeTeo.35Seo.65 where a small fluctuat- 
ing magnetic moment was found below T c .— This system 
again shows similar features in the Raman spectra as ob- 
served in FeSei^i 

Our most intriguing observation in FeSe is the on- 
set and development of a highly polarized low frequency 
background below 150 K which becomes gap-like under 
temperature lowering. This feature exists only in the 
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TABLE I. Results of phonon-mode calculations in FeSe and a comparison with data from Raman scattering experiments. 



Modes 


Experimental data 




PA/nmm at 300 K 




PA/nmm 


295 K 


7 K 


z Se = 0.2624 


z S e = 0.2493 


z S e = 0.2363 


setting 






not optimized 


optimized 


optimized 








M_Fe = 


Mp„ =2ub 


M_Fe = 


E 9 






124.3 


131.6 


143.0 


A lg (Se) 


179.8 


182.5 


178.3 


204.7 


210.6 


Bi 9 (Fe) 


193.9 


206.5 


218.8 


206.6 


230.8 


E u 






231.6 


255.3 


269.3 


E fl 






260.2 


294.7 


309.6 


A 2u 






276.6 


301.8 


293.4 



a&-polarization channel of scattering (B2 9 -channel in the 
tetragonal setting) . Meanwhile, we see an almost temper- 
ature independent background in all other polarizations 
(see Fig. 0}. To be sure with the sign and value of the re- 
spective Raman tensors we checked this response in RR 
and RL polarization geometry (not shown here) and no 
antisymmetric contribution into Raman tensors has been 
found. Thus we conclude that this background does not 
come from pure magnetic excitations or magnetic fluctu- 
ations. As seen in the Fig. [5]c the quasi-elastic part of 
the afe-Raman spectra starts to develop around 150 K and 
reaches a characteristic maximum of intensity at T s with 
a following decrease towards low temperatures. Note, 
that the order parameter rjs ~ u a {, (type of lattice dis- 
tortion) of the tetra-ortho phase transition has the same 
symmetry as the excitations seen in the afo-polarization 
channel of scattering. Due to this symmetry and charac- 
teristic increase at T s we ascribe the quasi-elastic part of 
scattering to the fluctuation of the u a b lattice distortions. 
Its existence in the Raman spectra in a wide temperature 
regime far above and far below of T s points to a very un- 
usual short-ranged lattice instability which persists even 
in the tetragonal phase. The rhombic instability of the 
tetragonal phase is usually considered a fingerprint of 
the electronic nematic fluctuations in underdoped iron 
HTSC42 

The gap-like feature shows a gradual increase of the 
frequency of its most intensive spectral weight starting 
from 150 K and up to 60 K. The frequency becomes sat- 
urated at lower temperatures at the energy ~ 450 cm -1 . 
A very similar ~ 290 cm energy gap with a similar- 
gap onset below 130 - 140 K has been observed recently 
in FeSe using time-resolved optical spectroscopy^ The 
gap size there was restored from indirect measurements 
and modeling which might explain the difference from 
our direct observation of the gap value. 

The gap onset accompanies a striking reduction in op- 
tical reflectivity and an increase of the Hall coefficient^ 
both signaling strong alterations in the electronic orbital 
and spin subsystems of FeSe below 150 K at the Fermi 
level. Furthermore, the hump excitations we observe 
have the symmetry of transitions between manifolds of 
electronic intermediate spin states and low and high spin 
states, that are allowed to be observed in ab polariza- 



tion. Thus, the gapped manifold of the intermediate spin 
states as well as accompanying orbital degrees of freedom 
do not participate in the formation of the paramagnetic 
ground state. This might explain the absence of orbital 
frustration in the FeSe paramagnetic ground state which 
we detect as a very narrow phonon linewidth. Interest- 
ingly, some alteration in the whole spin states manifold as 
well as in the appearance of the spin state gap have been 
supported also by observation of an anomalous change 
of the slope in the temperature dependence of the Se z 
coordinate that happens around 150 

Our findings also lead back to the question of the mag- 
netic origin of the structural tetra-ortho phase transition 
in the "11" family. In spite of its occurrence in FeSe with- 
out the onset of magnetic order like in FeTe, it occurs in 
FeSe on the background of some alteration in the mag- 
netic subsystems. Note, in previous studies the existence 
of a structural phase transition similar to FeTe in non- 
magnetic FeSe was the reason to decline the hypothesis 
about its magnetic origin^ 



V. SUMMARY 

Our experimental results as well as lattice dynamical 
calculations reveal a temperature evolution of the para- 
magnetic state in FeSe from the high iron spin state 
(S = 2) to the low iron spin state. We establish a high 
sensitivity of the Bi 9 (Fe) phonon mode, that is a charac- 
teristic mode of the FeSe4 tetrahedron layers, to the Fe 
spin state. Furthermore, we demonstrate the relevance 
of spin-orbital frustration despite a considerable metal- 
licity of these systems. Our scenario and the comparison 
of FeSe with FeTe underlines the complexity of Fe high 
temperature superconductors with respect to their lattice 
and electronic degrees of freedom. 
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